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Abstract 

We study the excitation spectrum and ground state properties of the two- 
leg S = ^ quantum spin ladder with staggered dimerization. Two massive 
phases, separated by a critical line are found, as predicted by previous analy- 
sis, based on the non- linear sigma model (NLSM). We have used dimer series 
expansions, exact diagonalization of small clusters and diagrammatic analysis 
of an effective interacting Bose gas Hamiltonian, obtained by using the bond 
operator representation for spins. We compute the excitation spectrum in 
the phase, dominated by strong rungs in the parameter regimes far away and 
close to the point of instability. The exact location of the phase boundary is 
determined. 
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I. INTRODUCTION 



There is a great deal of current interest in quasi one-dimensional quantum spin models, 
such as spin ladders, spin chains with dimerization and (or) frustration, and various gen- 
eralizations of the above. The studies of these models have been mostly triggered by the 
experimental discoveries of several spin-Peierls and spin-ladder compounds [|J. The theo- 
retical efforts have been focused mainly on the nature of the ground state, the excitation 
spectrum, as well as thermodynamic properties of these models. The system of coupled 
quantum S = | chains (spin-ladder) was found to behave differently depending on whether 
the number of chains is even or odd For an even number, the excitation spectrum is 
generically gaped and short-range correlations dominate in the ground state. The single 
5=1 quantum spin chain with dimerization exhibits quite similar properties, since it is 
known that dimerization produces an energy gap and destroys the quasi long-range order of 
the integrable uniform chain. 

In a recent paper, Ref. 0, M.A. Martin-Delgado, R. Shankar, and G. Sierra have pro- 
posed a model of coupled spin chains with staggered dimerization which possesses a rich 
phase diagram, quite different from that of the simple spin ladder (no dimerization). Con- 
sider the simplest version of the model, containing two coupled S = | Heisenberg chains 
(Fig. 1.): 

h = j±Y1 s m • S ^ + EE J ( 1 + s(-iy +a }s a , ■ s a>i+1 (i) 

i i a=l,2 

Here a = 1,2 is labeling the two chains (legs of the ladder). The coupling J_l > is the 
inter-chain (rung) interaction and J(l ± 5) > are the nearest neighbor interactions along 
the chains (all couplings are antiferromagnetic). Notice that the dimerization 5 is staggered 
in both directions (along and perpendicular to the chains), i.e. the dimerization in each 
chain is in antiphase with its neighbor chain. The model Eq.(l) was analyzed in Ref. H by 
mapping onto a non-linear sigma model (NLSM). In the NLSM approach the value of the 
topological angle 9 determines the nature of the possible phases of the model M. For 6 = 
the model is massive while it is massless for 6 = (odd multiple of n). For the model at hand 
two massive (gaped) phases were found in the parameter plane (S,J±/J) with a massless 
(critical) line separating them, determined by the equation: 

And 

Q = - — — = vr A5 = 2 + J L /J (critical line). (2) 

2 + jj_/2j 

A symmetric branch with 5 — > — 5 (8 = 3n) also exists but we will assume 5 > from now 
on. The two massive phases separated by this line are basically phases where the dominant 
correlations are along the rungs or the stronger bonds along the legs. The existence of a 
critical line is quite unusual since individually both interactions, J± and 6 are gap producing. 
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It is clear from Fig. 1. that there are two points (5, J±/J) = (0,0) and (1,2) where the 
system is a single uniform chain and hence integrable and massless. The NLSM approach 
Eq.(2) does not correctly reproduce the point (0, 0) which is quite natural since it is expected 
to work in the vicinity of the point (1, 2) only. It has been conjectured in Ref. [Q] essentially 
from continuity, that a critical line connects the two integrable points. 

Let us mention that upon introducing an additional spin alternation |4| and (or) consid- 
ering ladders of spin S = 1 or higher || the phase diagram of the system will be even more 
complex. The basic origin for this complexity however is still the staggered dimerization in 
the rung direction. Thus, form now on we will consider only the "minimal" model Eq.(l). 

The purpose of the present work is to analyze the excitation spectrum and the phase dia- 
gram of the model Eq.(l) at T = by a combination of numerical and analytical techniques. 
In Sec. II we present exact diagonalization results for the energy gap and the spin-spin cor- 
relations, confirming the existence of two gaped phases separated by a gapless line. In Sec. 
Ill we develop two kinds of strong-coupling approaches in the phase dominated by strong 
interchain (rung) correlations: (a) dimer series expansions, and (b) diagrammatic treatment 
of the effective Hamiltonian, written in terms of bond triplet operators. The spectrum near 
the critical line is calculated and the relevant interactions, responsible for the closing of the 
gap are identified. Sec. IV contains our conclusions and discussion of future prospects. 



II. EXACT DIAGONALIZATIONS 



First we present results from finite lattice diagonalizations of the Hamiltonian Eq.(l), 
for systems of up to N = 24 spins. Periodic boundary conditions are used throughout. For 
reasons of symmetry the number of spins must be a multiple of 4, and most results are 
obtained from extrapolation of 12, 14, 16, 20, 24 spin systems. We find empirically that the 
ground state is a singlet S to t = whereas the first excited state is a triplet S to t = 1. The 
computed singlet-triplet gaps are extrapolated to the thermodynamic limit N — > oo via the 
ansatz: 

A B 

AE N = AE^ + — + (3) 

N NlnN V ; 

where A and B are constants. The last term allows for the logarithmic corrections aris- 
ing through conformal invariance in spin chains and effectively accounts for the residual 
curvature observed in direct ratio plots. 

In Fig. 2. we show the variation of energy gaps versus the interchain coupling J±/J for 
various N, for a fixed value of the dimerization 5 (in this case 8 — |). It is clear qualitatively 
that the variation in AEn is small for both small and large values of J±/J whereas the 
variation is large near Jj_/J ~ 1.2 where the individual curves have a minimum. The dashed 
lines, based on the extrapolated formula (3) clearly show the gap vanishing at a single point 
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on the J±/J axis. Similar results are obtained for other values of 5 although for small 5 
the convergence is poorer. This is consistent with the increasing quantum fluctuations and 
larger finite size corrections expected near the uniform limit. This analysis leads to Fig. 3. 
which shows the estimated critical line in the (5, J±/J) plane. Comparing with Eq.(2) one 
can see that the NLSM result can not be trusted numerically, but the existence of a gapless 
line is correctly predicted. 

In Fig. 4. we show various spin correlations versus Jj_/J for 6 — |. As intuitively 
expected the phase with J±/J > 1.22 is characterized by strong inter-chain (rung) and 
weak intra-chain correlations. When J±/J decreases and ultimately crosses the critical 
point, the behavior of the correlation functions is reversed, signaling a transition into a 
phase, dominated by strong intra-chain correlations on the stronger bonds. We find that 
the spin correlations as well as the excitation gap change continuously through the quantum 
transition point. 

III. STRONG COUPLING EXPANSIONS 

In this section we present analysis based on strong-coupling expansions around the limit 
J_l ^> J. Two approaches are used: the dimer series expansion and the bond operator 
mapping onto an effective interacting Bose gas. The two approaches are similar in spirit 
and produce quite similar result even though the technical details are different. 

A. Dimer series expansions 

We start with results obtained by the linked-cluster dimer series expansion method. The 
linked-cluster expansion method has been previously reviewed in several articles |5|-0, and 
will not be repeated here. The basic idea is to divide the Hamiltonian into two parts: 

H = H + V (4) 

where 

Ho = J± Si^ ■ S2,tj 

i 

(5) 

V = JJ2 J2 t 1 + ^(-l) i+ 1S a ,,, ■ S a ,i+1 
i a=l,2 

By fixing values of 5, we can construct an expansion in J / J± by taking Hq as unperturbed 
Hamiltonian and V as a perturbation. The zeroth order V = J = approximation corre- 
sponds to isolated dimers (singlets) on the rungs. Each singlet can be excited into a triplet 
state and thus the excitation gap is J±. Finite J introduces interactions between the rungs 
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which modifies substantially the spectrum. In order to consider the physically interesting 
coupling region Jj_/ J ~ 1 the perturbation series has to be developed to high order and 
then extrapolated to the relevant point. We have carried out an expansion up to order 
(J/J_l) 23 for the ground state energy, to order (J/J±) 13 for antiferromagnetic susceptibility, 
and to order (J/ J±) 11 for the triplet excitation spectra for several values of 5. There are only 
12 graphs that contribute to the ground state energy and dispersion, and 14 graphs that 
contribute to the antiferromagnetic susceptibility. The series coefficients are not presented 
here, but are available upon request. Then we have used integrated differential approximants 
and Pade approximants to extrapolate the series. It is sometimes useful to compare the 
present results with the spectrum of a spin ladder without dimerization (5 = 0), which can 
be found in Ref. g (see Fig. 4.). 

In Fig. 5. we present the triplet excitation spectrum u(k) for 8 = |, obtained by the 
dimer series expansions. The triplet gap A = uj(it) decreases with the decrease of J±/J, 
vanishing (within error bar) at the critical value J±/J = 1.23. By comparing the spectra 
of Fig. 5. with the corresponding spectra for 5 = (see Fig. 4. of Ref. || for the value 
J±/J = 2) one can see that the renormalization due to finite 5 is the strongest at k = n 
while in the vicinity of the point k = the energy is almost unrenormalized. This is due to 
the fact that the first three leading orders in J/J± do not change uu(k = 0) (see also Sec. 
III.B). 

A more accurate determination of the critical line can be achieved by the Dlog Pade 
approximants to the antiferromagnetic susceptibility series, and the results are presented 
in Fig. 3. The phase boundary obtained within the dimer series approach is in excellent 
agreement with the exact diagonalization results in the region J±/J > 1. For small values 
of Jj_ the convergence of the dimer series becomes poorer (larger error bars), as the critical 
point occurs at much larger J/J±. In Fig. 6. several critical spectra, for parameters on 
the massless line, are plotted. Notice that for 5 = 1, Jj_ = 2 J when Eq.(l) reduces to a 
"snake", i.e. one- dimensional Heisenberg chain, the dimer series reproduces correctly the 
spinon dispersion, which is known from the exact Bethe ansatz solution of the problem [JTOj] . 
The latter is given by the formula (for a chain with exchange J): u(k)/J = (tt/2) sin( 2L =^). 
Indeed, the lower curve in Fig. 6. practically coincides with this formula. 

Using the long series for the ground state energy, we were able to obtain the most accurate 
estimates of the ground state energy, as it was done before for the normal two-chain ladder 



ljj . For example for 5 = 0.5 and J / J± = 0.8 which is near the critical line, the ground-state 



energy per site is estimated to be 

E /NJ ± = -0.52655(2) . (6) 

Beside the above inter-chain dimer expansions about the limit J± ^> J, we can also 
construct another dimer expansion, the intra-chain dimer expansion, by taking the stronger 
intra-chain bonds (the bonds couping J(l + <5)) as unperturbed Hamiltonian, and the rest of 
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the bonds as perturbation, and carry out dimer series expansions in (1 — 8)/ (1 + 8) for fixed 
values of J±/ J(l — 8). The series have been computed up to order [(1 — 8)/(l + 8)} 11 for the 
ground state energy, to order [(1 — 8)/(l + 8)] 6 for the antiferromagnetic susceptibility, and 
to order [(1 —8)/(l + 8)} 7 for the triplet excitation spectra for several values of J±/ J(l — 8). 
With this expansion, one can study the excitation spectrum and the ground state properties 
for the parameters located within the intra-chain dimer phase. The gapless critical line 
derived from this expansion agrees well with that by inter-chain dimer expansions, but with 
much less accuracy. Here we will not present any detailed results from this dimer expansion, 
however the series are available upon request. 



B. Bond operator representation 



In this section we describe a mapping of the model onto an interacting Bose gas, which 
can be achieved by using the bond operator representation for spins [[H|. Denote by |0)j the 
singlet state formed by two spins on rung i. Then, in the strong-coupling limit J±/J^> 1 the 
excitations are well described by triplets, created by the (bond) operators t ia , (triplet) 4 = 
t ia \0)i, where a = x,y,z are the three components of the triplet and the usual bosonic 



commutation relations hold [fL^] . The Hamiltonian Eq.(l) expressed in terms of the triplet 



operators is [13 



H = \ J Aitai + \ (4Wl + 44+1 + h-C.) 

i,a,/3 



where we have separated the part due to dimerization: 



^a^cw+i^j^yi + h.C. 



+ [i «-> % + 1]| 



(7) 



(8) 



and Hu will be defined below. For 8 = the Hamiltonian (7) coincides with that of the 
simple two-leg ladder [0,0 • In order for the triplet operators to represent only the physical 
spin triplet states, they have to satisfy the on-site hard-core constraint: 44 = 0- This 
restriction on the Hilbert space can be taken into account by introducing an infinite on-site 
repulsion between the triplets, as discussed in our previous work fI5| , |l6|| : 



Hu = U 2^44^**°*' U °°- 



(9) 



i,a/3 



Our goal is to develop a diagrammatic treatment of the interactions in the Hamiltonian 
(7) in order to understand which terms are most relevant to the renormalization of the 
spectrum. The quadratic part (first line in Eq.(7)) is diagonalized in momentum space by a 
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Bogoliubov transformation t ak — » Ukt a k + Vktl l _ k , which leads to the spectrum uj\ = A\ — B\, 
with A k = Jj_ + Jcos(fc), B k = Jcos(k). The Bogoliubov coefficients are defined as: 
u\,v\ = ±l/2 + A k /2uj k . 

The spectrum of the Hamiltonian Eq.(7) without the dimerization term H$ was investi- 
gated diagrammatically in Refs. |1|[1(J]. Here we only summarize the results. The on-site 
repulsion Hjj gives the dominant contribution to the spectrum renormalization, while the 
two-particle inter-site interaction with strength J/2 (first term in the second line of Eq.(7)) 
is a relatively minor effect ||17|| . Even though the on-site interaction is infinite, the scattering 
amplitude T a p^s(k, uj) = T(k,uj)(S ai 5/3s + SasSp-y) is finite (as physically expected), which 
can be seen by resumming the ladder series, shown in Fig. 7(a) and setting U — > oo. The 



scaterring amplitude and the corresponding self-energy (Fig. 7(b)) are []l"5"l , |l6 



[r^r^y; — — ±- — , E [/ (^) = 4y«;r(H^-4 (io) 

q H H q 

Here, and in all future equations we set, for convenience, the number of rungs to unity. As 
emphasized in Ref. fll5 , |l6l , Eq.(lO) represents the dominant contribution from Hu, provided 



the density of triplets N t = (t a>i t at i) = 3^2 q Vg is small. Thus Eq.(lO) should be viewed as 
the first term in an expansion in powers of N t . For the simple ladder (5 = 0) we find (after 
solving the corresponding Dyson equation, see below) N t « 0.1 for J_l = 2J, and N t « 0.25 
for J_l = J. Therefore the dilute gas approximation is expected to work quite well even 
for J±/J ~ 1. Naturally when J±/ J — > the above picture, based on strong rungs, breaks 
down. Instead, the excitations in this regime should be viewed as weakly bound spinons, 
thus making the local triplets an inadequate staring point. We have estimated that the bond 



operator formalism describes well the excitation spectrum for J±/J > 1 115,18 . 

Next we turn to the three-particle term H$ which represents the dimerization part of 
the Hamiltonian. To lowest (one-loop) order there are contributions to the normal and 
anomalous self-energies, drawn in Fig. 7(c) and Fig. 7(d), respectively. The corresponding 
formulas are: 

Ej(k, uj) = A{8Jf y C ^ + ^ k - <k k ) D ^ + -,k- g ,k) + ( u~v \ 

g UJ ~ ^9+- - ^k-q [UJ^-UJ j 

where the following definitions are used 

C(k, p, q) = $(fc, p)u k u p + $(-g, p)u k v p + $(fc, -q)v k u p , (13) 

D(k,p, q) = —C(k,p, q){u <-> v}, $>(k,p) = k ~ smp). (14) 
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The function <&(fc,p) originates from the Fourier transform of the three-particle vertex Eq.(8). 
Equations (11) and (12) can be derived by evaluating all possible internal loops in the 
diagrams in Figs. 7(c), (d), where, as usual, lines with a single arrow stand for the normal 
Green's function G N (k,t) = —i(Tt a ^(t)V ak (0)) and lines with double arrows represent the 
anomalous Green's functions G A (k,t) = — i(TV a _ k (t)t' a fe (0)). Finally, the coupled Dyson's 
equations for the normal and anomalous Green's functions have to be solved leading to the 



result (more details can be found in [IS 



C N (h . .\ = U + A k + E N (k, ~UJ) 

l ' ' [u + A k + E N (k,-u)][u-A k -J: N (k,u)} + [B k + E A (k,u)}^ 1 ] 

where the total normal and anomalous self-energies are defined as 

T, N (k, u) = Z v {k, u) + (k, u), S A (fc, u) = Zj(k, u). (16) 

In order to get a feeling for the effect of H& on the spectrum, it is instructive to examine 
the second order perturbation theory result. It can be obtained by noticing that to leading 
order oj k w J±,u k « l,t> fc « — (J/2J±) cosk. Thus in evaluating the leading order correction 
arising from the dimerization term one can set v k = which means that only the normal 
self-energy contributes: 

5u k = Ef(k,u; k ) = 4(5 J) 2 V ^ + "' k - q) = ^(cos k - 1). (17) 

q ~ J± J± 

One can see that to lowest order the renormalization of the spectrum for k = ir is the 
strongest while k = is not renormalized. For a fixed value of 5, the gap cj(7t) decreases 
with decreasing J±/J and vanishes at a critical ratio, signaling an instability of the rung 
dimer phase. 

We have found the spectrum numerically by solving for the poles of the exact Green's 
function Eq.(15), thus resumming the perturbation series to all orders. The results for 8 = ~ 
are presented in Fig. 5. The diagrammatic results are in good agreement with the dimer 
series expansions, however with decreasing J±/J the disagreement increases. We would like 
to point out that even though the two methods have the same starting point (strong rungs) 
and typically agree quite well in the region of strong coupling, for smaller coupling the dimer 
series expansion performs better. The reason is that the series expansion takes into account 
all terms to a particular (finite) order, while the diagrammatic approach sums up only the 
most dominant subclass of diagrams. The small parameter N t controlling the diagrammatic 
expansion grows with decreasing J±/J and at J±_ = 1.23 J (the critical point, as estimated 
by the dimer series and exact diagonalization) has the value N t = 0.3. This relatively large 
value of the triplet density is related to the considerably decreased magnitude of the spin- 
spin correlation function on the rungs |19|] (see the (1,2) correlations in Fig. 4.). Therefore 



the agreement between the diagrammatic analysis and the dimer series is, in fact, even better 
than expected. 
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IV. SUMMARY AND DISCUSSION 



In conclusion, we have investigated the excitation spectrum and the ground state phase 
diagram of the quantum spin ladder with staggered dimerization and have found that a 
transition takes place form a phase dominated by strong rung correlations into a phase, 
characterized by intrachain dimers. The two phases are separated by a gapless line. The 
location of the boundary, the triplet excitations, as well as the spin-spin correlations in the 
rung dimer phase were calculated by dimer series expansions, diagrammatic analysis of the 
effective hard-core Bose gas of rung triplets, and exact diagonalizations. 

We have identified the three-particle scattering H$, in terms of the rung triplets, to be 
responsible for the instability of the rung dimer phase (closing of the gap). It is interest- 
ing to note that this is indeed the term which, in the continuous field theory formulation, 
contributes to the topological angle 9 and leads to the transition. Indeed, from the repre- 
sentation of the spin operator on a given rung i in terms of the triplets |R| one can easily 



see that Si — S 2 ~ t, + tj ~ Si + S 2 ~ ztj x t; ~ lj. These are precisely the two vector 
fields used in the derivation of the non-linear sigma model 0. Therefore the terms of the 
form cf> i • lj + i which appear in our lattice Hamiltonian H$ would be precisely the ones which 
give a non-trivial 9 after the continuum limit is taken. 

There are several additional aspects of the problem that can be addressed within the 
strong-coupling formalism described in this paper. As suggested in Ref. the spinons 
become unconfined on the gapless line separating the two phases of the model (see Fig. 3.). 
Since in the two massive phases confinement definitely takes place, we would expect that 
the triplet spectral weight decreases as the critical line is approached and ultimately goes to 
zero. In addition, at the point of deconfinement the same energy is required to create triplet 
and singlet excitations. Therefore, close to the critical line, a low-energy singlet excitation 
must appear in the spectrum. Within the strong-coupling approach a singlet can be viewed 
as a collective bound state of two triplets |15| . While in the spin ladder without dimerization 
the singlet bound state is high in energy (small binding energy), in the present model we 
expect it to become a truly low-energy state (strong triplet binding). We leave these issues 



for a future study |20| 



While this manuscript was being prepared for publication we became aware of two recent 
preprints devoted to the same model ^J. Both of these works have studied the phase 
diagram numerically via exact diagonalization and their results for the location of the critical 
line are in agreement with ours. 
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FIGURES 

FIG. 1. The spin ladder with staggered dimerization. 

FIG. 2. Energy gap for different system sizes (solid lines) and its extrapolation to the thermo- 
dynamic limit (dashed line). 

FIG. 3. Phase diagram of the model. Open circles are the exact diagonalization data, where 
the circle diameter represents the error in the determination of the critical point. Solid squares 
with error bars are the critical points obtained by Dlog Pade approximation of the dimer series. 

FIG. 4. Spin-spin correlations (SfSj) (where (i,j) are defined in Fig. 1.) obtained by exact 
diagonalization of a system of N = 20 spins. The dashed line represents the location of the critical 
point. 

FIG. 5. Triplet energy spectrum obtained by the dimer series expansion (solid squares con- 
nected by solid lines) for J±/J = 2,1.43,1.23 (upper, lower and middle curve at k = it, respec- 
tively). The corresponding spectra, obtained from the poles of Eq.(15) are drawn, respectively, 
with solid, long-dashed and short-dashed lines. 

FIG. 6. Same as Fig.5. for several values of (S, J±) on the gapless line. 

FIG. 7. (a) Resummation of the ladder series for the scaterring amplitude T(k,cu), where k(u>) 
is the total incoming momentum (energy), k = k\ + fe. (b) The normal self-energy corresponding 
to T. (c) Lowest order diagrams for the normal self-energy arising from the three-particle term H$. 
(d) Same as (c) for the anomalous self-energy. The dots represent all other possible diagrams with 
different internal loops. 
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